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We construct a non-perturbative approach based on quantum averaging combined with resonant 
transformations to detect the resonances of a given Hamihonian and to treat them. This approach, 
that generalizes the rotating-wave approximation, takes into account the resonances at low field and 
also at high field (non- linear resonances). This allows to derive effective Ifamiltonians that contain 
the qualitative features of the spectrum, i.e. crossings and avoided crossings, as a function of the 
coupling constant. At a second stage the precision of the spectrum can be improved quantitatively by 
standard perturbative methods like contact transformations. We illustrate this method to determine 
the spectrum of a two-level atom interacting with a single mode of a quantized field. 

PACS numbers: 03.65.-w, 02.30.Mv, 42.50.Hz, 42.50.Ct 



I. INTRODUCTION 

Some important features of classical and quantum sys- 
tems are determined by resonances of the system which 
can not be treated by perturbative approaches. In the 
vicinity of resonances the perturbative formulas display 
small denominators that lead to the divergence of the per- 
turbative expansions. A widely used model that incor- 
porates a one-photon resonance is the Jaynes-Cummings 
Hamiltonian extracted from the full dressed Hamiltonian 
that describes a two-level system coupled with a single 
mode of a quantized field Q . Its counterpart for an inter- 
action with a semi-classical laser field is the RWA Hamil- 
tonian (rotating-wave approximation) 

In this article we give a systematic method that al- 
lows to construct effective Hamiltonians and determine 
their spectrum by treating the resonances with an adap- 
tation of resonant transformations that were introduced 
in Ref. Q in the context of laser-driven quantum sys- 
tems in the Floquet representation. The goal is to ob- 
tain the spectrum for a whole interval of values of a pa- 
rameter like the coupling constant. This is needed e.g. 
in applications where the coupling changes adiabatically 
corresponding e.g. to envelopes of laser pulses or to 
transversal spatial profiles of cavity fields. The method 
is based on the detection of resonances by a projector de- 
rived from Quantum Averaging (QA). We illustrate it on 
the problem of a two-level atom interacting with a quan- 
tized field and show that a treatment of all the relevant 
resonances of the system in a given range of parameters 
allows to reproduce with good accuracy the spectrum of 
this system. The treatment of the resonances yields the 
qualitative structure of the spectrum - the crossings and 
avoided crossings - as a function of the coupling constant. 
Once this main structure is obtained, one can systemati- 
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cally improve the quantitative accuracy of the spectrum 
by applying perturbative methods. We use contact trans- 
formations with a Kolmogorov-Arnold-Moser (KAM) it- 
eration j^, that are particularly efficient due to its su- 
perconvergent properties. 

The paper is structured as follows. In Section (nj we 
describe the method of resonance analysis and the con- 
struction of effective Hamiltonians. Section ITTll contains 
the presentation of the model and some preliminary con- 
siderations. In Section Hvl taking into account the res- 
onances of this model in the weak-coupling regime, we 
extract the effective Hamiltonians by quantum averaging 
techniques and resonant transformations. In the weak- 
coupling regime we have to iterate this procedure several 
times to derive the essential structure of the spectrum 
in larger ranges of the coupling constant. In Sec|Vl we 
extract the effective Hamiltonians in the strong-coupling 
regime where the qualitative properties of the spectrum 
can be globally obtained by some preliminary unitary 
transformations and one resonant transformation which 
treats the zero-field resonances. We obtain an accurate 
approximation valid for all values of the coupling con- 
stant that contains all the qualitative structure. Finally, 
in Sec I VII we give some conclusions. 



II. PRINCIPLE OF THE METHOD 

We consider a Hamiltonian H = Hq + eV where Hq 
is the unperturbed Hamiltonian, eV is the perturbation 
and e is an ordering parameter. The first analysis of this 
problem is in terms of perturbation theory: we look for a 
KAM-type unitary transformation e*^*^ close to the iden- 
tity that allows to reduce the order of the perturbation 
from e to e^: 

e-'^He'^ = Ho + eD + e^V2. (1) 

eD is a remaining term of order e that satisfies [Hq, D] = 
0. The unknown W and D are solutions of the following 



2 



equations [1,13 



[Ho,W] + V = D, 
[Ho,D] = 0. 



(2a) 
(2b) 



The remaining perturbation of order is given by 

^'^2 = E S ((™ - l)L"^-'V + L^-'D) , (3) 



m=2 



where is defined as 



LwB = [B, W] . 



(4) 



The solutions of Eqs.© can be written in terms of aver- 
aging 



D = V = nH„V:= Um - f dse-'"°'Ve'"°'' 



(5a) 
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where z/ labels the different eigenvalues e'^^'^ of ffo, and j 
is a degeneracy index which distinguishes different basis 
vectors of the degeneracy eigenspace. The operator 
H-Hq is the projector on the kernel of the application A ^ 
[Hq^A\. We remark that the integral representation of 
Z?, W in Eqs. lO can be also well-defined in cases where 
Hq has a continuum spectrum. The units are chosen such 
that h = 1. In the following discussion, we do not write 
explicitly the ordering parameter e. 

A resonance is defined as a degeneracy of an eigenvalue 
E^u^ of Hq and is said to be active if the perturbation V 
has nonzero matrix elements in the degeneracy subspace 
{v:j\y\v:j') for some j,/. Otherwise the 
resonance is called passive or mute. An active resonance 
renders W arbitrarily large close to the degeneracy and 
makes the perturbative expansion diverge. The method 
we present here, is a construction designed to avoid such 
divergences. We remark that the concept of resonance is 
defined intrinsically for _ffo, while the distinction between 
active and passive depends on the relation between 
and V . The analysis of the resonances involves thus three 
aspects: 

• Decomposition of the Hamiltonian into H — Hq + 
V. Different decompositions can be considered for 
different regimes of the parameters of H . 

• Determination of degenerate eigenvalues of Hq. 

• Detection of the resonant terms in the perturbation 
V that couple these degenerate eigenstates. 



The resonant terms of V can be detected by projectors of 
type Uhq that extract a block-diagonal part of V relative 
to Hq, where the blocks are generated by the degeneracy 
subspaces. In absence of active resonances, when all the 
eigenvalues of Hq are non-degenerate or when the reso- 
nances are mute, the matrix representation of TIhoV is 
in fact diagonal in the eigenbasis of Hq. In presence of 
active resonances, the block-diagonal effective Hamilto- 
nian that takes into account the considered resonance of 
the original Hamiltonian can be written as 



H 



off 



Ho+IIh„V. 



(6) 



We will call the transformation that diagonalizcs H'^^ 
Resonant Transformation (RT). The Hamiltonian H = 
i/cff^ (^V-HhoV) is transformed under RT (denoted 7^) 
as follows: 



Hi=n'^Hn = 7^^i^°^7^ + 7^t(y-HH„v^)7^ 

Vu (7) 



=: 



where Hf'^ is defined as the new renormalized reference 
Hamiltonian and Vi is the new perturbation. If H'f^ + Vi 
does not have any other active resonance in the consid- 
ered range of the coupling constant, we can at a second 
stage improve the spectrum by a KAM-type perturba- 
tive expansion which is expected to converge. If there are 
other active resonances, we have to iterate the renormal- 
ization procedure by applying another RT. We remark 
that there are cases of multi-photon resonances where the 
resonant terms appear only after applying one or several 
contact transformations. 



III. DESCRIPTION OF THE MODEL AND 
PRELIMINARY CONSIDERATIONS 



We consider as an illustration a two-level atom inter- 
acting with a single mode of a quantized field described 
by 

H = u}{a^a+l/2)iS)t2+'^t'»a^+g{a + a'^)<^c7x, (8) 
where a, are the annihilation and creation operators 



for the field mode with the commutation relation I 



1 = X]^ol"')("^l ' '^zjO'x are Pauli matrices and I2 is 
the 2x2 identity matrix. Here oj is the frequency of field 
mode, luq is the energy difference of the two atomic states 
and g is the dipole-coupling between the field mode and 
the atom. This Hamiltonian acts on the Hilbert space 
K. = J-®TL where = is the Hilbert space of the atom 
generated by |±) (eigenvectors of CTz) and T is the Fock 
space of the field mode generated by the orthonormal 
basis {\n) ; n = 0, 1, 2, • • • }, n being the photon number 
of the field. 

For this system there is a parity operator 



n=0 



(-l)'»(n|®a. 



(9) 
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with the properties 

[P,H]^Q, P^P\ 1^ = 1® 12. (10) 

As a consequence, the eigenstates of H can be separated 
into two symmetry classes, even or odd, under P: 

P|<^„,±) = ±\^n.±), i?|^„,±) - £;„,±|0„,±). (11) 

The parity operator also commutes with any operator 
that depends only on = a^a and az- 

In spite of the simple form of ((HJ, its exact solutions 
are not known. This can be related to the fact that the 
classical limit of this model is non-integrable 0- This 
model is ofgreat interest as a physical model in quantum 
optics IE 0| and quantum chaos 0, hj. Some 
approximate solutions of this model have been studied 
among many others in 14, 15j using different formalisms. 

The conceptual framework for the solution of this sys- 
tem based on the construction of unitary transformations 
can be described as follows: First, we decompose the 
Hamiltonian in two terms as H = Hq + V. Depending on 
the considered ranges of the parameters of the system, 
different decompositions may be considered. Hq is a pri- 
ori an operator that is a regular function exclusively of 
the operators iV and cr^ . The operators N and can be 
considered in the present model as quantum analogues 
of classical global actions 16], and Hq can be labelled 
integrable. The perturbation V contains functions that 
involve also the other operators a, a^, CTj,, ctj^. The goal is 
to determine a unitary transformation U, that should be 
expressed in terms of well-behaved regular functions of 
a, a'' ,ax,a'y,<Jz, such that: 

t/t {Ho{N, a,) + y(a, a\a,,ay,az)) U - H'{N, a,), 

(12) 

where H' is a regular function / exclusively of the ac- 
tion operators iV,cr^: H'{N,az) = f{N,c7z). With this 
transformation the eigenvectors of H can be expressed 
as \4>n,±) = U{\n) (g) |±)) and the corresponding eigen- 
values as En ± = f{n,±l) where N\n) = n\n) and 
a,|±)=±|±)'. 

We remark that in our context the important property 
for singling out the operators N, is that they com- 
mute with each other and their spectrum and eigenvec- 
tors are explicitly available. The question of whether for 
a given model there exists a regular unitary transforma- 
tion U that accomplishes the above requirement is, to 
our knowledge, an open problem. 

Most of the perturbative approaches can be interpreted 
as methods to find approximations of the transformation 
U . The presence of resonances is one of the central diffi- 
culties in the construction of f7, as will be made precise 
below. In this paper we discuss an iterative approach 
that consists of constructing first some approximations 
of U that take into account the dominating effects of a 
certain number of resonances. The transformations in- 
volved in this stage are far from the identity and have 
a clearly non-perturbative character. Once we have a 



transformation that takes into account the main effect of 
a set of resonances that are relevant in a considered inter- 
val of the coupling constant 5, a perturbative approach 
(like the KAM, Van Vleck, or other types of the contact 
transformation) can be applied to improve the approxi- 
mation quantitatively. The transformations involved in 
this second stage can be considered as deformations of 
the identity, since they can be written in the form . 
This stage cannot be implemented if the resonances are 
not taken care of beforehand. Indeed the perturbative 
formulations diverge close to resonances due to the ap- 
pearance of small denominators as can be seen in Eq. 
(Hb). 

As in classical mechanics, the construction of the trans- 
formation U leading to a Hamiltonian that contains only 
action variables can often be considered in two steps: 
U = UiU2- In the first step, that is called reduction, the 
Hamiltonian is transformed by C/i into a form that con- 
tains functions of CzjCx, cry and N , but not of a and a^. 
The degree of freedom of the field is made trivial and the 
number of non-trivial degrees of freedom is thus reduced 
by one. When we apply this reduction to the effective 
Hamiltonian lO, we obtain a reduced effective Hamilto- 
nian. We remark that in the literature, this "reduced 
effective Hamiltonian" is often called simply "effective 
Hamiltonian" . In the second step, the reduced Hamilto- 
nian is transformed under U2 into a form that contains 
functions of only A^ and az- For the model JSJ, the re- 
duction step corresponds to diagonalization in the Fock 
space and the second step corresponds to diagonalization 
in the atomic Hilbert space which in this case is trivial. 
The construction of the RT is based on this reduction 
procedure. 

IV. EFFECTIVE HAMILTONIANS IN THE 
WEAK-COUPLING REGIME 

In this section we consider the Hamiltonian © at res- 
onance ojQ = Lo 'm the weak coupling regime, so that H 
can be decomposed as follows: 

H = Ho + V, 
HoiN,az) = t^(Ar+l/2)®l2-H yl^fJ^, 
V{a,a\ax,g) = .g(a + a^) (g) CTj,. (13) 
The eigenvalues and eigenvectors of Ho are: 

4% = uj{n + l/2)±uJo/2, 
|<^i%) ^ K±) = |n)®|±), 

For loq = Lo there is a one photon resonance which corre- 
sponds to the degeneracies = -E^+i^.. The degen- 
eracy eigenspaces are spanned by the vectors and 



4 



\4>n+i -)■ The resonant part of V is obtained by jSJ-a): where 



Vres-^'^HaV = ^ ( +) (n, + 1 V |n + 1 , -) (n + 1 , - 1 
n=0 

+ |n+l,-)(n+l,-|y|n,+)(n,+|) 
where we have used the relations 

OO CXj 

V" + 1 |?i)(7i + l|, at = V Vn + l + 



E 

ri=0 



n=0 



(16) 

The effective Hamiltonian containing the one-photon res- 
onance is the so-called Jaynes-Cummings Hamiltonian 
that can be written as 



ifo = Hjc = Ho + nHoV = Lu{N+ 1/2) ®t2 



-1 ®(Jz+g 



at 



(17) 



Hjc is a good approximation of ((HJ for low energies in 
the limit g <^ lJ{)^ \ lo ~ lJ{)\ loq. In this limit, the 

so-called comiter-rotating terms 5 ^ ^ ^ '^^^ 

carded (rotating- wave approximation). H can thus be 
written as H = HQ^{N,a,a\ax,(Jy]g) + {V — UffgV). 
Next we transform Hq^ by a resonant transformation TZi 
to a regular function of exclusively the action operators 
N, ■ Every resonant transformation is performed in 
two steps. To diagonalize in the Fock space (the 
reduction step of the RT denoted Ri) we define a trans- 
formation in such a way that the following condition is 
satisfied: 



(18) 



where / is a regular function of N which has to be de- 
termined. We require furthermore that r\HqRi stays a 
function of only N and cr^. A suitable transformation 
satisfying these conditions is 

^ /(aat)-V2« 0\ /Er=ol">(" + l| 

y 1/V 1 

This transformation is not unitary but isometric 17]: 
?t_,.. ptp_,._f I0)(0| 



RiR^ — Ijc , R^Ri — Ifc — 



/ ' 



(19) 
1: 

(20) 



where we have used the identity at(A^+l) ~ 1— | 
Applying this transformation on the resonant term gives 

R\VresRi = g a^aa^y^/'^a ®ax= gVN (g) (21) 
and H is transformed under i?i as 

HR,^R\HRi^cuN(g)l2+gVN<»<7x+g(^''^ ^' 

(22) 



A = a(aat)-i/2a = ^V^TTT |n)(n -I- 2|, (23) 



with the properties: 

AA^ = aat 

To each eigenvector 
of since: 



^t A = a^a- a-f |0)(0|. (24) 
of H corresponds an eigenvector 



HR,R\\(j)) = R\HRiRl\(f)) = \r\ 



(25) 



We remark that R\\(j)) ^ V|(^) G /C. Every eigen- 
value of the original Hamiltonian H is also an eigenvalue 
of the transformed Hamiltonian Hr-^ . However since 
i?i|0, = 0, there is a difference in the spectrum be- 
tween H and Hr^ : Hr^ has an extra zero eigenvalue with 
eigenvector |0, The spurious eigenvalue can be de- 
tected and eliminated after applying the transformation. 
Indeed, since |0, is not coupled to any vector in its or- 
thogonal complement, one can eliminate it from the rest 
of the calculation by taking the projection of Hr^ into the 
orthogonal complement i^J^l,J_(o,+) = fj_(o,+)-H^fii^J-(o,+) 
with P_i(o,+) = Ik; — |0, +)(0, +|.This difference between 
unitary and isometric transformations was not taken into 
account in 18] in diagonalizing the Jaynes-Cummings 
Hamiltonian. 

The second step of the RT is the diagonalization of 
rIhI^^Ri = ojN (g)l2 + VN (g) (T^ in the atomic Hilbert 
space. This can be performed by a 7r/2 rotation around 
the y-axis: 



with the properties 



(26) 



(27) 



However, since the spurious eigenvector |0, can 
be separated and |0, — ) is already an eigenvector of 
R\Hq^Ri, the transformation T must be applied only on 
the subspace with n > 1 photons. The complete trans- 
formation (denoted Ti) reads thus 



where 



Po = |0)(0|, P^o = 5]|n)(r 

Tl = l 



(28) 



(29) 



Applying Ti gives 

M 

h[°'^ (N, a,;g) + Viia,a\a,,ax,ay;g), (30) 



Hi := tIr\HRiTi 
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FIG. 1: Comparison of exact numerical eigenvalues (dashed 
lines) of © for one-photon resonance cj = luq with the ap- 
proximate ones (solid lines) obtained after (a) 1 one-photon 
RT given by (I34L (b) 1 one- photon RT plus 1 iteration of 
KAM-type perturbative expansion. The divergence observed 
around g/uiQ = 0.3 in panel (b) is due to the active nonlinear 
resonances of occurred at the degeneracies marked by 

circles in panel (a). One can see clearly that the locations of 
these resonances depend on n according to Eq. I|85|l . 



with 



(0) 



A±o - A{^ - A^o - A{, 



V, = 9/2 ( ' " 



_9_[ |2)(G| \ 

^[\0){2\ -|2)(0|-|0)(2|;' 



where 



Aio = ^'±o^^±o-^\A^|?i)(« + 2|, (32) 



and use has been made of the relations 

APo = PoA^ = 0, PaAP^o = |0)(2|. (33) 

The first RT, is thus the combination of i?iTi. Since the 
transformation i?i dresses the upper atomic state by (-1) 
photon 0, TZi = RiTi can be called a one-photon RT. 

H^^^ is in fact the diagonalized Jaynes-Cummings 
Hamiltonian in the resonant case with the eigenvalues 



e: 



tL±)i9)^^n±gV^, n = 0,l,2,.... (34) 



(0) 



The eigenvalues and therefore the degeneracies of Hi 
depend on the coupling constant g. For small enough 



g and low energies, H^'' does not have other degenera- 
cies besides the ones at 17 = for which the new per- 
turbation Vi does not have resonant terms, and we can 
apply KAM-type transformations to improve quantita- 
tively the precision of the spectrum by iteration. A sin- 
gle KAM transformation (which is essentially equivalent 
to second order perturbation theory) already gives quite 
good precision, as shown in Fig. I^b) for g/ujo < 0.25 
for energies smaller than lOwo. If we take large enough 
g or larger energies, we encounter new resonances which 
appear at some specific finite values of g. These reso- 
nances are called field-induced resonances or nonlinear 
resonances. For larger values of the coupling (g/uJo ~ 0.3 
for the shown energy interval in Fig. (Q-b)), where 
we encounter nonlinear resonances, the KAM iteration 
diverges. The eigenvalues of h[^^ are degenerate at 



g„ - ^/(V^+ V^TTT) as E^°l^+-^{gn) = i?i;'(Ui,_)(5n)- 
But the corresponding resonant terms in Vi are zero due 
to parity (mute resonances). The next degeneracies ap- 
pear at 



(0) 



5„ = 2uj/{V^+Vn + 2), 



as 



E 



(0) 

!,(",+) 



E. 



(0) 

l,(n+2,-) 



i9n), 



(35) 



(36) 



which have been marked by circles in figure l^a). All 
the other resonances are mute. There is an infinite fam- 
ily of nonlinear resonances located at different values of 
the coupling gn . We observe from H35|) that for higher en- 
ergies the nonlinear resonances appear for arbitrary small 
coupling (lim„__>oo <7n = 0). We can extract the resonant 
terms corresponding to the whole family in a single step 
by working with the combined projector J2n^H'°\g )■ 
The resonant terms in Vi corresponding to the degenera- 
cies are 

V - 9 f A^o 

~ ^{0 |2)(0| + |0)(2|)' (3^) 
and the new effective Hamiltonian is thus 

Hf ^LuN(g)l2+ gVN^cr, + n^(«)(g„)^i- (38) 

n 

To diagonalize H^^, it can be decomposed according to 
three orthogonal subspaces: 

Hf = P(o,2,-)i^r%,2,-)+^^(o,+)i^r%,+)+^'±i?f^± 



= HfPio,2,-) + HfP^o..+) + HfP± 



(39) 



where the projectors, which commute with Hf, are de- 
fined by 
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(0,2, 



Pi 





|0)(0| + |2)(2|y' 



1/C - -P(0,2,-) - -P(0,+) 



(0, 



|0)(0| 




En=i\n){n\ 




which leads to 



H^Pio. 



HfP± 



0, HfPio.2,-) 


gfO ^±0 



(2^-.gV2)|2)(2|-4(|2)(0| + |0)(2|) 



V2 
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E,T=i,„#2 l^^>(^ 



Hi^P{o,2,-) can be directly diagonalized by 





P(0,2,-) - P(0,2,-) I Q 



COS 



?(|2)(2|-|O)(O|)-sin0(|2)(O| + |O)(2|) 



(0,2,-), 



where the angle 9 is defined by the relation 



tan 29 = 



gV2 



2w -5V2' 



0<.<I, 



and the corresponding eigenvalues are 



E' 



icff 

1,(0,+) 
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The reduction step of the second RT to diagonalize Hi^P± in the Fock space can be defined as 



R2,± ■■= P± 
with the properties 



Er=i,#2l")(^ 



R2,±rI^ = p±, rI^R2,± = p±- 



|1)(1| + |2)(2| 




(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



Equation (|45|l shows that i?2,_L dresses the upper atomic state by (-2) photons. Therefore 7?.2,± can be called a 
two-photon RT. Since i?2,j.|l, +) = = i?2,±|2, +), the spectrum of rI j^Hf^P±R2,± has two extra zero eigenvalues 
relative to the spectrum of Hf^P^. Applying i?2.± gives 



^2,-L-^i*^^-^^2,± 



Er=3("-2)|n)(n| ^ 

Er=i,#2"ki>("l 



Er=3 - 2 \n){n\ 








■T,n=i.,^2Vn \n){r 



5/2^ Vti^^ \n){n\ (g) ax 



(47) 



n— 3 



Combining the transformations on the different subspaces we can write the transformation that diagonalizes Hf^ in 
the Fock space as 



i?2 — R2,± + P(0,2,-) + P(o,^ 



(48) 



At the right hand side of 14711 . the three matrices have entries that commute with each other so we can diagonalize 
the sum of them in the atomic Hilbert space (the second step of Tl2.± ) as if they had scalar entries. The eigenvalues 
of Rl^HfP±R2.± are thus: 



0.5 1 0.5 1 0.5 1 




0.5 1 0.5 1 0.5 1 

g/«o 9/®0 

FIG. 2: Comparison of the exact numerical eigenvalues (dashed lines) of ||HJ for one-photon resonance to — uo with the 
approximate ones (solid lines) obtained respectively after (a) 1 one-photon RT given by (I34I I. (b) 1 one-photon RT plus 1 
two-photon RT given by l|49|l . (c) 1 one-photon RT plus 2 two-photon RT, (d) 1 one-photon RT plus 3 two-photon RT, (e) 1 
one-photon RT plus 4 two-photon RT, (f) 1 one-photon RT plus 4 two-photon RT plus 1 iteration of KAM-type perturbative 
expansion. The divergence of the KAM transformation observed close to g/iJo = 1 in panel (f) is due to the presence of active 
resonances at larger values of g. 



pcff 

-'^l,(n>3,±) 



E' 



l,(n=l,+) ~ 

■9, 



pcff 



= 0, 



1/2 



(49) 



As it can be seen from ()49f) there are two extra zero eigen- 
values which have been added by i?2.i to the spectrum 
of Hf. 

Figs. l|21-a,b) compare respectively the exact spec- 
trum of H calculated numerically with the spectrum of 
Hf = Hjc given by ^ and of Hf given by gni),131- 
The crossings of the exact spectrum are all among the 
eigenvalues with different parities. It is found that the 



spectrum of Hf coincides with the exact one only in the 
range of quite small coupling. The spectrum of Hf has 
been modified with respect to the one of Hf by trans- 
forming the encircled crossings between eigenvalues with 
the same parity into avoided crossings in the small g re- 
gion. This procedure to treat resonances can be iterated 
to take into account other resonances appearing at larger 
values of 5. Figs. |21-a,b,c,d,e) show how the combination 



8 



of a one-photon RT and consecutive two-photon RTs hft 
the artificial degeneracies (marked by circles) of the effec- 
tive Hamiltonians. The successive steps which we have 
implemented numerically, transform eigenvalue crossings 
into avoided crossings. We observe that these RTs also 
produce an improvement of the approximations of the 
spectrum. Fig. |21-f ) shows the effect of a KAM transfor- 
mation after the fourth two-photon RT which improves 
quantitatively the result of Fig. (|2J-e). The divergence of 
the KAM transformation close to g = 1 in Fig. ||21-e) is 
due to the presence of active resonances at larger values 
of g. 



V. EFFECTIVE HAMILTONIANS IN THE 
STRONG-COUPLING REGIME 

In this section we use quantum averaging techniques 
and RT to obtain the effective Hamiltonians of |(HJ) in the 
strong-coupling regime. We derive a formula that repro- 
duces the spectrum quite accurately in the whole range 
of g and for all energies. We consider the Hamiltonian 
(|5]) in the strong coupling regime (7 ^ cjq > , which 
suggests to decompose the Hamiltonian H as 



H 
Ho 
V 



Ho + V, 

cj(iV+ l/2)(» I2 
1 a,- 



9[a 



^1, 
2 



(50) 



that can be interpreted as the system of a quantized field 
plus the coupling term perturbed by the two-level atom. 
We remark that in this decomposition, Hq contains all 
the unbounded operators of the complete model and that 
the perturbation y is a bounded operator. In this case 
Hq{N, a, flt, az, Cx! g) is integrable since we can explicitly 
transform it into a form involving a regular function ex- 
clusively of the action operators N, Oz (given below in 
Eq. H55|l ). To transform iJo to a function of action oper- 
ators, first we diagonahze the term g{a -|- a^) (g) cr^ in the 
atomic Hilbert space by the transformation (|26l) : 



Next we apply a second unitary transformation 



T'^UT = w(iV+l/2)®l2+.g(aW)«)cr^ — ^l^cr^. (51) 



U = 











(52) 



to transform uj(N + 1/2) (g) I2 + g(a + a^) ® Oz into a 
function of only N ^ Oz (in this case only of N) : 



Hi := U^T'^HTU 





-2i(at-a) 







(53) 



where use has been made of the commutation relations 
among a, a\ N and the Hausdorff formula: 



^^Ce-^=C+[i?,C] + ^[i3,[S,q] 



1 



We decompose Hi as 
Hi = Hf^ +Vi, 

Vi = ~uJa/2 



2!' 



c^(iV + l/2)-^ 

UJ 

J 



(54) 



(55) 



The effective Hamiltonian of the system for strong- 
coupling regime can thus be written as 



iff -ijf^+n^(o,Fi 

r(0) 



(56) 



The eigenvalues of H{ have a two-fold degeneracy for 
every value of n as 



<J,,±)=-(n-fl/2)-^ 



The average of Vi relative to H^^^ is thus 

00 

n<)'^'i = ^{|n,+)(n,+|yi|n,-)(n,- 

n=0 

+ \n,-){n,-\Vi\n,+){n,+\} 

OC 



(57) 



(58) 



with 

In 



[me 



+ 2s 




\n) 




^/ 




V {n-i)\ 



E 



(+2g/^) 



(n — i)! 



(j!)2 {n~jy. 



E 

3=0 



(59) 



where the Ln are the Laguerre polynomials. We remark 
that in the limit of large photon number (n 00), /„ 

can be expressed as a zero-order Bessel function Jo(^^j^) 
[lol |. Hi can be reorganized as 

Hi = Hf+iVi-U^mVi), 
Hf = ('c^(iV + l/2)-^') 12- yi^^a,. 



(V n V) - ^0 G-F 
(yi-H^(o,yi) _ ^ 



(60) 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



which is the same result obtained in [Hill Hi by other 
methods. Figure Q compares the exact numerical spec- 
trum of JHl with the approximation H65|l for the resonant 
case cj = cjo- One can see that for large enough g, the 
formula (|65|l reproduces well the spectrum. It is not very 
accurate for small values of g because of the presence of 
the one-photon zero-field resonances that we analyze as 
follows. In the limit 17 — > 0, we have 

V'r"".(V-a.)°0°')- 

Thus degeneracies occur as 



FIG. 3: Comparison of exact numerical eigenvalues (dashed ^(0),s^0 _ ^(0),g^0 (Q^) 

lines) of as a function of the coupling constant in the 2,(n,+) ~ 2,(ri-l,-)' v ' 
resonant case {uj = ojo), with the approximate eigenvalues 
(solid lines) obtained from 1651 . 



They are made active by the resonant terms of V2 



9-0. 



where 



G = e^('^'-'^), F = Y,fn\n){n\. (61) 



n=0 



HI can easily be diagonalized by applying the trans- 
formation H2()|) that diagonalizes cTj.: 



(62) 



with 



Hf^ = T^HfT 



u}{N + 1/2) - ^ ) ® I2 - ® (T„(63) 
UJ J 2 



and 



V2 = T\Vi-n^mVi)T 

= -^o/^[ _G + Gt -G-Gt + 2Fj-(64) 
The eigenvalues of iJj"'' are therefore 
41 ±) - -(-+1/2) -'-T ^e-^V.^i„( V ). (65) 



^2, res ^ 2 



= -.9 



at 
a 



(68) 



The transformation (the reduction step of the RT) which 
transforms this resonant term to a regular function of N 
is 



Ri 



1 \ _ f ^ 
(aat)-V2„; - 1^0 E,T=ol")(" + liy ' 

(69) 

with the properties 



RiR\ = lie, R\Ri — lie 





|o)(oiy' 



(70) 



We remark that the definition of Ri depends on the 
type of resonant terms. The reduction step of the RT 
presented here is different from (|19|l . The Hamiltonian 
transformed under this RT has an extra zero eigenvalue 
corresponding to spurious eigenvector |0,— ), while for 
the Hamiltonian transformed under (|19|l . the extra zero 
eigenvalue corresponds to |0,-|-). Applying Ri on H2 
gives 
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FIG. 4; Comparison of exact numerical spectrum of l|HJ (dashed lines) as a function of the coupling constant in the resonant 
case {lu = luq), with the quite accurate result 17511 which has treated the zero-field resonances by a RT (solid lines). 



H3 RIH2R1 = (ujN 



Next, we take i/g"' = ujN 1^ I2 and the rest of H3 as V3. Since H^^'' has a two- fold degeneracy as E'. 
the average of V3 relative to H^^'' is thus 



I2 + RIV2R1 



f (i-E,T=o/nl«)("l) 

(0) 



(71) 



E. 



(0) 



3,(n,+) -^3,(n,-)' 



-^e-^^'/-'L^^\C^)\n){n\ -(| + f^)(l - |0)(0|) + f ^^=1 fn-i\n){r 
where we have used the relation [Til 



with "''(2;) the associated Laguerre polynomials. The new effective Hamiltonian can thus be written as 



H. 



(0) vz- 



(72) 



(73) 



(74) 



Since all the entries of commute with iV, it can be diagonalized in the atomic Hilbert space as if its entries were 
scalars. The eigenvalues of H^^ are thus 



El 



eff 

3,(0, 



^3,(0,+) - 2 " ~ ~ 2 ' " 



-^3,(n>l,±) 



o2 

u} 4 ' 



nto re I Ln[^r) - Ln-l[^r) 



± i 
2 



n 



1/2 



(75) 



The zero eigenvalue is the extra spurious one that has pares the exact numerical spectrum of ||SJ) and the ap- 
been added by the RT to the spectrum. Fig. Q com- 
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proximation H75() which has treated the zero-field reso- 
nances by a RT. The figure shows that treating all the 
active resonances of the system allows to obtain all the 
qualitative features of the spectrum in the whole range 
of the coupling constant and for all energies. At a second 
stage, since we have treated all the active resonances, we 
can improve further this spectrum quantitatively by a 
KAM-type perturbative iteration. 

VI. CONCLUSIONS 

We have presented a non-perturbative method based 
on the quantum averaging technique to determine the 
spectral properties of systems containing resonances. It 
consists in the construction of unitary or isometric trans- 
formations that leads to an effective reduced Hamilto- 
nian. These transformations are composed of two qual- 
itatively distinct stages. The first one consists of non- 
perturbative transformations (RTs) that are adapted to 
the structure of the resonances. Their role is to construct 
a first effective Hamiltonian that contains the main qual- 
itative features of the spectrum - crossings and avoided 
crossings - in a given range of the coupling parame- 
ter. The diagonalized form of this effective Hamiltonian, 
which depends parametrically on the coupling constant, 
is then taken as a new reference Hamiltonian around 
which one can apply perturbative techniques to improve 
the quantitative accuracy of the spectrum. We formulate 
the perturbative approach in terms of a KAM-type iter- 



ation of contact transformations. Similar results can be 
obtained with other formulations of perturbation theory. 

We have illustrated the method with a model of a two 
level atom interacting with a single mode of a quantized 
field. The method can be applied to more general systems 
with several field modes. It can also be adapted to the 
treatment of semi-classical models in which the field is 
described as a time-dependent function. 

We have analyzed the resonances in two regimes of 
weak and strong coupling. The results we obtained in the 
weak-coupling regime can be expected to be applicable to 
quite general models. The analysis of the strong-coupling 
regime of this model leads to results that are valid for all 
values of the coupling and for all energies. The possibil- 
ity to obtain such a global result is due to a particular 
property of the model, and one cannot expect to obtain 
it for general models. The particular property is that 
the part we selected as the reference Hamiltonian Hq in 
the strong-coupling regime contains all the unbounded 
operators of the complete model and is explicitly solv- 
able. The term that was left to be treated by RT and 
perturbation theory is a bounded operator. 
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